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Abstract 

Hamiltonian lattice QCD with Wilson fermions is investigated systematically by strong- 
coupling expansion up to the second order. The effective Hamiltonian is diagonalized 
by Bogoliubov transformation. The vacuum energy, chiral condensate, pseudo-scalar and 
vector meson masses are calculated. The comparison with the unitary transformation 
method by Luo and Chen is also made. The method discussed in this paper has potential 
application to QCD at finite density. 
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1 Introduction 



Lattice Gauge Theory (LGT) proposed by Wilson |J is the most promising non-perturbative 
method for QCD. In comparison to other techniques, the advantage of LGT is that there is 
no free parameter when the continuum limit is taken. There are two different approaches: 
Lagrangian formulation and Hamiltonian formulation. Lagrangian formulation is conve- 
nient for Monte Carlo simulations. Hamiltonian formulation is useful for analytical calcu- 
lations H H, £|, H H, [7j . Although the standard lattice Lagrangian Monte Carlo method has 
been successful in many aspects of LGT, it unfortunately breaks down at finite chemical 
potential due to the "complex action problem" . On the other hand, lattice QCD at finite 
chemical potential formulated in the Hamiltonian approach does not encounter the com- 
plex action problem. Recently, Gregory, Guo, Kroger and Luo || investigated the finite 
density Hamiltonian lattice QCD with naive fermions, using the unitary transformation, 
variational method and Bogoliubov transformation method developed by Luo and Chen 



It is well known that naive fermions experience "species doubling". Kogut-Susskind 
fermions and Wilson fermions are the most popular ways towards solving of the prob- 
lem. Although the doubling problem has not yet been completely solved, Kogut-Susskind 
fermions are frequently used for exploring the spontaneous breaking of the U(l) chiral 
symmetry. For Wilson fermions, the doubling modes are removed, but chiral symmetry is 
explicitly broken and fine tuning of the fermion mass remains to be done; This method is 
very popular in spectrum computations, due to the existence of the flavor symmetry. 

In this paper, we study the Hamiltonian lattice QCD with Wilson fermions by strong 
coupling expansion method. The effective Hamiltonian is Fierz-rearranged, re-expressed 
by meson operators, and then diagonalized by Bogoliubov transformation. The vacuum 
energy, chiral condensate, pseudo-scalar and vector meson masses are computed. 

The rest of the paper is organized as follows. In Sect. 0, we obtain the effective 
Hamiltonian using strong coupling expansion. The physical results are presented in Sect. 
^| and Sect. [|. The conclusion is provided in Sect. [S[ 

2 Effective Hamiltonian at strong coupling 

We begin with the (d+l)-dimensional lattice QCD Hamiltonian with Wilson fermions 
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where m, a, r and g are respectively the bare fermion mass, lattice spacing, Wilson param- 
eter, and bare coupling constant. U(x, k) is the gauge link variable at site x and direction 
k, and if) is the four-component spinor. The convention = — 7^ is used. E a (x) is the 
color-electric field and summation over a = 1, 2,..., 8 is implied. U p is the product of gauge 
link variable around an elementary plaquette and in the continuum, it represents the color 
magnetic interactions. However, at the strong coupling, the color magnetic energy (the last 
term) can be ignored. 

We want to diagonalize H so that the fermion field ip can be expressed in terms of up 
and down 2-spinors £ and rf 

*M = ( X) ) ■ (2) 

The bare vacuum state |0) is defined by £|0) = 77 10) = E?(x)\0) = 0. Since the up and 
down components are coupled via the 7^ matrices, the bare vacuum is not an eigenstate of 
H. 

For the convenience of strong coupling expansion, we introduce the un-perturbative 
term W and perturbative term W, i.e., 

^H = W + W, (3) 

g 2 

where 

a g x ti 

W = \Y,^( x )(~ r + ^)U(x,k)ij(x + k). (4) 

g x,k 

The ground state energy is g 2 e/2a, where for 1/g 2 << 1, e can be approximated by 

e=(0\^H\0} = e^+e^+e^. (5) 
9 

e^°\ e^ 1 ' and are the zeroth, first and second order corrections: 

e {0) = Wb|0>, 
£ « = (Q\W\0), 

From (||) and (^), we obtain 
e« = 0, 
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(0|^^(x)V(^)|0) = -2N e N f N,—. (7) 
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Table 1: T matrices and coefficients for H*ff . 

Here iV s , iV c and Nf are respectively the number of lattice sites, colors and flavors. M 

vd 

M = m + —, (8) 
a 

and l/(2Ma) is an analog of the k parameter in Lagrangian LGT. Using 

W o W\0) = \W ,W]\0) + WW \0) 

= \W o ,W]\0)+eV>W\0), (9) 

the second order correction becomes 

e {2) = (0|W ^~^y |0) = -^(0|^(1 - |0}(0|) W|0> = -i-(0|W|0>, (10) 

where C N = Tr J2 a X a X a /N c = (JVg - 1)/(2N C ), i.e. the Casimir invariant. From <g) and 
(fTo|) , we have 

^ = ^ (0) -(0|^|0) = (0|^ e// |0). (11) 

Therefore, the effective Hamiltonian is 

Heff = Mj2Hx)tP(x) + Hff + H e k ff . (12) 

X 

This equation will serve as our starting point, where 

Kr 2 

H r ff = _ onrE^./i( I )^./i( a; + ^OVWafa + k)ip Clj2 (x), (13) 

and 

K 

H k = 7r^Y,^c 1 jA x )ik^c 2 jA x + k )^c 2 jA x + k )ik^ ci j 2 {x). (14) 

laN c X)k 

K = 1/ ((7 2 Cat) is the effective coupling of the four-fermion interactions. The color index c 
and flavor index / are explicitly written, and summation over repeated index is implied. It 
is easily proven that the contribution from terms proportional to 0(r) vanishes. 
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Table 2: T matrices and coefficients for H%. . 



In order to express the effective Hamiltonian by meson operators, we perform a Fierz 
transformation so that fll3|) and ([14]) are rearranged as, 



and 



o-e// 



Kr 2 d Kr 2 

£ ^(x)V(x) + — £ L A V& (x)V>} 2 (x + fc)^ (x + fc), (15) 



K 



iaN, 



Y,L A ^ h {x)Y A ^ h {x)^ h {x + fc)rA^A(x + k). (16) 



c x ,fc 



where the matrices Ta and their coefficients L' A and are given in Tab. [1] and Tab. [2], 
and summation over the index A is understood. Note our basis of the sixteen Ta matrices 
is the same as the standard reference fLl| , but it is different from SmitQ. 
Substituting flTSP and (fEB) into flI2|), we have 

K 

a 



H, 



eff 



+ 



K 



laN r 



^ ^ (r 2 + l)V' t (x)V'(a;)V' t (x + fc)^(x + fc) 



'c a; fe=±j 

+ (r 2 — l) , 0'''(a;)74'0(a;)'0 (a; + k)^7p(x + k) 

- (r 2 — l)-0T( a ;)7 5 '0(a;)'0''(a; + k)j5ip(x + k) 

+ (r 2 + l)^ t (x)7 4 75-?/>(£)'?/> t ( a; + fc)7 '47 ' 5 i/j (x + fe) 

+ (r 2 + (1 - 25 fcj -)) ^ t (a;)747iV'(x)^ t (x + k)^Jjip(x + fc) 

- (r 2 - (1 - 24j)) ^(x)^(x)^(x + k)jjip{x + k) 

- {{I - 25 k j) + r 2 ) il)\x)-i i ajil)(x)il)\x + fc)7 4 ay0(;E + fc) 

+ ((1 - 26 k j) - r 2 ) ip i (x)a j ip(x)^(x + k)ajip(x + k) 



(17) 



where o,j = £jj 1 j 2 ljiJj2- This effective Hamiltonian describes the nearest-neighbor four 
fermion interactions. For r — * 0, H e ff is the same as that from the unitary transformation 
method using by Luo and Chen|J. 

Define the pseudo-scalar meson operator II and the vector meson operator V as [Bf 



U(x. 



hh 
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x 



U^(x) f2 fl = 2 ^— ^ 2 (x)(l +74)75^/1^), 

Then Eq. (|T7|) can be re-expressed as 
H eff = (—2N c NfV) (m + —^- \ + 2M £ (l I ; , /: I, ,,(,•) • Vj^ 

V a I xjxfro 

2Kr 2 d 

— E (njyiWAW+^W^A^ 

- — £ nt i/2 (x)n /2/l (x + fc) + n /l/2 (x)n} 2/i (a; + A;) 



+ V kh WVjhh (% + k) + V jhh {x)V] hfl {x + k) 
K 
2a 



+ I; E (nU^) 11 /*/^ + fc ) + n /i/ 2 (^)n/ 2/l (x + k) 



x,k 

+ (^A/2 + k ) + W*)%2a(* + *)) (1 - 2**j) 

V/1/2 as/i/aj 

= (-2N c N f N s ) (m+^J+H ti + H v , (19) 
where the contributions of the pseudo-scalar mesons are vector mesons are respectively 

tfn = (-™£* + m) E n} 2/l (.)n /l/2 (,) 

\ a ) xhh 

Kr 2 



E (I"/: /,( '')'"/:-/:(•'• + *) + Uf:fM)U M: (r + fc 



and 



^ a xfifok 

+ ^ E (4 l/2 (x)4 2/i (x + fc) + n /l/2 (x)n /2/l (x + fc) 

+ E 4 2/i (a;)n /l/2 (x), (20) 

^ = 2M-^— + — £ 4 2/l (x)% l/2 (x) 
\ ' ' s/i/2i 

6 



+ ^ E (4 1 / 2 (^)4 2 / 1 ^ + fc ) + ^ 1 / 2 (^)% 2 / 1 (^ + ^))(i-2^)- 

xflhkj 

3 Meson masses 

After a Fourier transformation 

p 

Hu in (^) becomes 

#n = f^(l-r 2 ) + 2M) J] 4 i/2 (p) a/2/l (p) 

\ / p/l/2 



Kr 



2 



£ E ( a /!/2(P) a /2/l(P) + a (P)/l/2 a /2/l(P)) E COS Pi a 



/l/ 2 P 



+ — E (4i/ 2 (P) a / 2 /i(~P) + a /i/ 2 (-p) a / 2 /i(p)J E cos Pi a - 

The Bogoliubov transformation 

a(p) — > a(p) coshw(p) + a^(— p) sinhw(p), 

cr(p) — > a*(p) cosh-u(p) + a(— p) sinhw(p), 

o(— p) — > a(— p) cosh -u(p) + cr(p) sinhw(p), 

a ' ( — p) — > a ' ( — p) cosh u (p) + a (p) sinh u(p) , 

will lead to the diagonalization of ifn, if the parameter w p satisfies 

tanh2w(p) = — — — cos p q a, 
G i q =i 

where 

G - * 
a 



Kd, 2N Kr 2 * 

M H ( 1 — r ) > cos p„a 

a a ^ 
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Therefore, after the Bogoliubov transformation, the diagonalized Hn becomes 



H n = Gi i 1 -tanh 2 2w(p)J 2 a) ih (p)a hh (p) 
p/1/2 



1 — tanh 2 2u(p) 



2G 2 r 2 * 

— — 2^ cos Pi a 



9=1 



(27) 



From (^3), (pq) and (p7[), we obtain the dispersion law for the pseudo-scalar mesons, i.e. 
the relation between the energy En of the pseudo-scalar meson and momentum p: 



E u = Gx (l - tanh 2 2u(p)) 5 = G 1 



,2God 



Gi 



M + 1 -r 2 



cosp ? a 



9=1 



x 



M + — (1 - r 2 ) - — (r 2 + 1) cosp q a 



9=1 



(28) 



For naive fermions (r = 0) in 3+1 dimensions, according to fl28|) , we reproduce the result of 



Luo and Chen||: E n \ Pq=0 
the dispersion law tells us 



m n oc m, i.e. the PCAC theorem; In the chiral limit m — > 0, 



4i^ 2 



9 - (^cosp 5 a) 
9=1 



(29) 



One sees that there exist two zeroes at p— (0, 0, 0) and p = (tt, it, it). This is a well known 
"doubling problem". For Wilson fermions (r ^ 0), the doubler modes are removed, but the 
chiral symmetry is explicitly broken. In order to define the chiral limit, one has to fine-tune 
M — > M c so that the pseudo-scalars become massless. From (|28|), we get 



Mr 



6Kr 2 



(30) 



According to this formulae and (0), near the chiral limit, the mass of a pseudo-scalar 
behaves as 



-^nU=o 



2AK 

4(M - M c f + (M - M t 



24 K 
a 



(31) 



which is the PCAC relation for a pseudo-scalar in the Wilson fermion case. 
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We now consider Hy in ([Hp similarly. After a Fourier transformation 

p 

and a Bogoliubov transformation 

bj (p) — > 6j (p) cosh f j (p) + £>t ( — p) sinh ■Uj (p) , 
&j (p) — *■ (p) cos h u j (p) + ( ~P) s i nn u j (p) j 



where 



we obtain 



tanh2u,(p) 



2G 2 
G 1 



cos p q a — 2 cospj-a 



£ G 1 (l -tanh 2 2^(p))H] /i/2 (p)6 j/2/l (p) 
pj/1/2 

Gi 



PJ 



1 - (l -tanh 2 2^(p) 



2G 2 2 
1™ 



^ COSp q a 



9=1 



Therefore, using (|30|) , the dispersion relation for the vector mesons is 



2K 

a 



3(1 + r 2 ) + (1 — r 2 ) cosp g a — 2 cos p^a 

9=1 



3(1 + r 2 ) — (1 + r 2 ) cosp g a + 2 cos p^-a 

9=1 



For naive fermions (r = 0) in 3+1 dimensions, from ([5t]) and fl3"6]). 

4i^ 2 



9 — ( cos p g a — 2 cospj-a 



(32) 



(33) 
(34) 



(35) 



(36) 



(37) 



where the chiral limit has been taken. One sees that there exist six zeroes on the boundary 
of the first Brillouin zone, p = (it, 0, 0), (0, it, 0), (0, 0, n), (it, it, 0), (0, it, it), (tt, 0, n), due to 
the "doubling problem". For Wilson fermions (r 7^ 0), the doubling problem is avoided, 
and the mass for a vector meson is 



my 



Ej\ v= $ = d (1 - tanh 2vj(p)) | p=0 



G l 



2G 2 
~G~i 



(d-2) 



\p=0i 



(38) 
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where 



has been used. In the chiral limit M = M r 



4v/ V 

m v = K, 

a 

which is independent of the Wilson parameter r. 



(39) 



4 Vacuum energy and chiral condensate 

The vacuum energy is the vacuum expectation value of the Hamiltonian 

= (n\H\n) = (o\H eff \o) 



-2N c N f N s ) [ M + 



d 



g 2 C N a 
1 - (l - tanh 2 2mQo) ) 2 ) -\ 



a 



q=l 



Pj 



1 — (l — tanh 2 2vj(py) 2 ) + -p^—r 2 cosp q i 



2G 2 2 , d 



q=l 



In the r = limit, this agrees with Luo and Chen||. 
The fermion condensate 

(W) = -2N c NfN s + (^)n + 
can also be computed using the Feynman-Hellmann theorem. Here 

"ggi l £(Gl-(2G 2 j:Ucosp q a) 2 ) dG 
dm 2 • !G\ - (2G 2 Ef =1 cos^a) 2 



N 2 
f \ ~* 

2 ^ 

v 



(40) 



(41) 



9m 



2 2r 2 ^ cosp g a 



iV 2 iV s /•*/«» . 



(2vr) d J-Va 



AT? AT /-tt/o 



if V^a 



L M+ ^ (1 _ r2) _i^ E ^ cosp9a 



(27r) d y-Tr/a 



'i-[ 



~ 2^=1 cos P(3 a 



W£i£ + Jfd(i_ r 2)_2fi£y-" iCO sp„a J 

a a v J a Z— / g— 1 ^ y 



A^iV, fir/a 



(27T) d y-Tr/a 



d d P 



'1- 



6r 2 +d(l — r 2 ) — r 2 ^ _ 1 cosp g a 
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W)1 



N]N S 



it /a 
-7r I a 



d d P 



1 



^[J2 q = l COS Pl a - 2cOS P3 a 



\ l M+^(l-r»)-^E; =1 



cos p 9 a 



(27T) 



E 



7r/a 
— 7r/a 



1 



r X)g = i cos Pg a ~ 2cos Pj a ]2 
6r 2 +d(l— r 2 )— r 2 YJ _j cosp g a 



(42) 



Due to the Wilson term, ipip mixes with the identity operator []T0 

/.continuum „J,„i,lattice _|_ (jl J 



(43) 



This relation is useful for comparing the lattice result with the continuum theory when 
1/g 2 » 1, but the computation of the mixing coefficient C 1 is quite involved. 



5 Discussions and Expectation 

In the preceding sections, we have investigated (d+l)-dimensional Hamiltonian LGT with 
Wilson fermions. Using a strong-coupling perturbative expansion, the effective Hamiltonian 
in the strong-coupling regime has been obtained and diagonalized exactly by Bogoliubov 
transformation. Some interesting physical results for the vacuum energy, meson masses 
and fermion condensate have been obtained. For r = 0, our results reduce to those of Luo 
and Chen|5|]. We will apply these techniques to lattice QCD at finite chemical potential. 
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